A New Mechanism for the a to u Martensitic Transformation in Pure Titanium 
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We propose a new direct mechanism for the pressure driven a—*Lu martensitic transformation in 
pure titanium. A systematic algorithm enumerates all possible mechanisms whose energy barriers 
are evaluated. A new, homogeneous mechanism emerges with a barrier at least four times lower 
than other mechanisms. This mechanism remains favorable in a simple nucleation model. 

PACS numbers: 81.30.Kf, 64.70.Kb, 5.70.Fh 



Martensitic transformations are abundant in nature 
and are commonly used in engineering technologies 0. 
Materials from steel to shape-memory alloys are governed 
by their underlying martensitic transformations . The 
pressure driven a (hep) — > u> (hexagonal) transformation 
in pure titanium, discussed here and reviewed extensively 
by Sikka et al. |3j, has significant technological implica- 
tions in the aerospace industry because the u phase for- 
mation lowers toughness and ductility. This transforma- 
tion was first observed by Jamieson and has since 
been extensively studied using static high-pressure 
and Shockwave methods . Because of experimental dif- 
ficulties in directly observing martensitic transformation 
mechanisms, they are usually inferred from the orienta- 
tion relationships between the initial and final phases. 
Such an approach may result in multiple transformation 
mechanisms for any given set of orientation relations, 
and requires one to guess the appropriate transformation 
mechanism. Thus, despite several attempts, the mecha- 
nism for this transformation is still unclear. 

We calculate the energy barrier for homogeneous trans- 
formation for different titanium a— >lu transformation 
mechanisms and compare the values using a simplified 
nucleation model. We systematically generate and sort 
possible a— >lu mechanisms by their energy barriers. A 
new direct mechanism emerges whose barrier is lowest 
both homogeneously and when considered in a simple 
nucleation model. 

Figure n shows our new low energy barrier mechanism 
for the a—>uj transformation in Ti, called TAO-1, for "Ti- 
tanium Alpha to Omega." This direct 6-atom transfor- 
mation requires no intermediary phase, and has small 
shuffles and strains. The 6 atoms divide into a group of 
4 atoms that shuffle by 0.63 A and 2 atoms that shuf- 
fle by 0.42 A. Combining these shuffles with strains of 
e x = 0.91, E y = 1.12, and e z = 0.98 produces a final ui cell 
from our a cell. The original a matrix is then oriented 
relative to the lo matrix such that (0001) Q || (0111)^ and 
[1120] Q II [Olllju,. These orientation relations are seen in 



some experiments, but not others 

Our mechanism identification method matches possi- 
ble supercells of a and lo to determine the lattice strain, 
and atom positions to determine the necessary internal 
relaxations; similar to While there are infinitely 

many possible supercells, we consider only 6 and 12 
atom supercells with principal strains less than 1.333 and 
greater than 0.75 = 1.333 . For each supercell, there 
are multiple ways to shuffle the atom positions from a 
to lo; we consider mechanisms where, relative to the cen- 
ter of mass, no atom moves more than 2.0 A. We also 
limit the closest distance of approach between any two 
atoms during the transformation. Each supercell has one 
mechanism with the largest closest approach distance; we 
consider all mechanisms for that supercell within 0.1 A 
of that value. 

To efficiently compare the enthalpy barriers of the pos- 
sible mechanisms, we use three methods of increasing 
computational sophistication and accuracy. 1. For each 
mechanism, we calculate an approximate barrier based 
only on the strain, using the elastic approximation. Such 
a calculation of the "elastic barrier" relies only on the 
elastic constants for each phase, and the supercell geom- 
etry [^J. 2. For mechanisms with a low elastic barrier, 
we calculate an energy landscape using a strain variable 
(e) and a single concerted shuffle variable (5). The en- 
ergies E(e , S) are calculated using a tight-binding (TB) 
model [lOj on a grid between a at E(0, 0) and u at 
£7(1,1). The energy at the transition state in this re- 
duced space gives a "landscape barrier." 3. For mech- 
anisms with the lowest landscape barriers, we calculate 
the exact barrier using the nudged elastic band (NEB) 
method ^lj together with our TB model or ab initio 
VASP simulations or both. During the NEB calcu- 
lation, the cell shape and size are allowed to change in re- 
sponse to the applied pressure. In general, we expect the 
elastic barrier to underestimate the TB landscape bar- 
rier, and the TB landscape barrier to overestimate the 
TB-NEB barrier. The most probable mechanism should 
have the smallest ab initio NEB enthalpy barrier. 
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(a) hep cell (Q) and shuffled pos (Q) (b) hep cell with omega bonds (c) shuffled hep cell (d) omega cell (after shearing hep) 




FIG. 1: Our proposed a—>u> transformation mechanism (TAO-1). The direct 6-atom transformation is visualized (i) as a 
two-step process, (ii) with 21 additional atoms outside the heavy gray parallelpiped supercell. In a, atoms 1,3,5 and 2,4,6 are 
in the A and B stacking planes, respectively; whereas in uj, atoms 1-4 are in the Wyckoff d position of space group P6/mmm, 
and atoms 5-6 are in the Wyckoff a position (a) The grey atoms in the a cell shuffle to new positions (white atoms), with 
atoms 1-4 shuffling 0.63 A, and atoms 5-6 shuffling 0.42 A. (b) The a cell is redrawn with the "bonds" of the uj structure. 
Step 1: Shuffling the grey atoms to the white positions, the a cell (b) produces a strained ui cell (c), contained in the same 
supercell. Step 2: Straining the supercell (c) e x = 0.91, e y = 1.12, and e z = 0.98 produces the final w supercell (d). The 
orientation relations connecting the a and ui supercells are: (0001) a || (0111)^ and [1120] Q || [0111]^. 



We calculate total energies using a TB model and care- 
fully converged ab initio calculations. The TB calcula- 
tions are performed with OHMMS and use Mehl and 
Papaconstantopoulos's functional form |l(J with param- 
eters that are refit to reproduce LAPW total energies for 
hep, bee, fee, omega, and sc to within 0.5 meV/atom. 
For each mechanism, we use a k-point mesh equivalant 
tol2xl2x8in hep with a 63 meV thermal smearing 
to give an energy con verg ence of 1 meV/atom. For our 
ab initio calculations |2jJ we include 3p electrons in the 
valence band, and use a plane-wave kinetic-energy cutoff 
of 400 eV and a 7 x 7 x 7 k-point mesh to ensure en- 
ergy convergence to within 1 meV/atom. We relax the 
atomic positions and the unit cell shape and volume until 
the total electronic energy changes by less than 1 meV. 

The elastic and landscape energy barrier calculations 
reduce our initial set of 977 candidate mechanisms to 
the 3 lowest barrier mechanisms. From the initial list of 
134 supercells, we reject all but 57 because their elastic 
barriers are large (greater than 100 meV/atom). For the 
remaining 57 supercells, we generate 977 mechanisms (6 
6-atom mechanisms and 971 12-atom mechanisms), and 
compute TB landscape barriers. Figure [21 shows that the 
clastic barrier underestimates the landscape barrier. We 
select the three mechanisms with the smallest landscape 
barriers — TAO-1, TAO-2, and Silcock — for careful study. 

Figure |31 illustrates the Silcock mechanism 0] . Her 
mechanism, determined from observed orientation rela- 
tions, involves significant atomic shuffle, relatively small 
strains, and is a direct transformation mechanism with 
no intermediate state. In each a stacking plane, 3 out 
of every 6 atoms shuffle by 0.74 A along [1120] Q , while 
the other 3 shuffle in the opposite direction [1120] a . This 
shuffle is accompanied by a strain e x = 1.05 along [1100] a 
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FIG. 2: Approximate elastic energy barrier (meV/atom) ver- 
sus landscape energy barrier (meV/atom) for 359 mechanisms 
whose landscape barrier is below 200 meV/atom. Of the 134 
supercells generated, 57 had low (less than 100 meV/atom) 
elastic barriers. From those 57 supercells, 977 mechanisms 
were generated; for those 977 mechanisms, we calculate a 
shear/shuffle energy landscape using our tight-binding param- 
eters. Of those 977, 359 mechanisms had an elastic barrier 
below 200 meV/atom. We select the three mechanisms in the 
bottom left corner — TAO-1, TAO-2, and Silcock — and calcu- 
late their true barrier using the ab initio NEB method. 



and e y = 0.95 along [1120] Q to produce a hexagonal ui 
cell with the correct c/a ratio. Despite the lack of direct 
conclusive evidence to support the Silcock mechanism, it 
is still widely invoked to describe the a—>-uj transforma- 
tion because it is a direct mechanism [l5|. 

The remaining two mechanisms, TAO-1 and TAO-2, 
are related to Usikov and Zilbershtein's proposed 
composite mechanism involving the thermodynamically 
unstable /3(bcc) intermediary phase, and proceeding as 
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FIG. 3: Silcock mechanism for a to ui transformation (single 
a stacking plane shown), (a) In each stacking plane, 3 out of 
every 6 atoms shuffle by 0.74 A along [1120] Q , while the other 
3 shuffle in the opposite direction [Il20] a . (b) This creates the 
(1120) W plane from (0001) Q . The a and d Wyckoff positions 
of PQ/mmm are also shown. This shuffle is accompanied by 
a strain e x = 1.05 along [llOO]^ and e y = 0.95 along [1120] a 
to produce a hexagonal ui cell with the correct c/a ratio. The 
orientation relations connecting the a and ui supercells are: 
(0001)a || (1120). and [1120] Q || [0001] u . 



TABLE I: Comparison of TAO-1, TAO-2, and Silcock mech- 
anisms. Energy barriers: Four different methods for calculat- 
ing the energy barrier for the three mechanisms are shown, 
from least accurate to most accurate. The elastic barrier only 
accounts for the strain in each mechanism. The landscape 
barrier uses a simple combined shuffle for each, and a tight- 
binding total energy. Finally, the NEB calculation is done 
with the tight-binding method and ab initio to accurately de- 
termine the barrier. Orientation Relations: The relative ori- 
entation of a to ui is shown for each mechanism. N.B.: Silcock 
and TAO-2 mechanisms have the same orientation relations. 
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a— +/3— >u). Using the ol-uj orientation relations, they pro- 
posed an a— >0 transformation via an inverse Burgers 
mechanism | l6j followed by a (3— >uj transformation via 
the collapse of 2 out of 3 (111)/? planes. This produced 
two unique mechanisms, called Variant I and Variant II, 
depending on the direction of {111}^ planes to collapse. 
The TAO-1 and TAO-2 mechanisms could have been con- 
structed from Variant I and II, respectively, by using 
the j3 intermediate system to construct a 6 atom super- 
cell. Allowing Variant I and II to relax away from the (3 
phase would have resulted in direct transformation mech- 
anisms. 

Table [I] summarizes the energy barriers for the three 
mechanisms of interest and identifies TAO-1 as having 
the lowest energy barrier. The TAO-1 mechanism pro- 
vides the best trade off between shuffle and strain; dur- 
ing the transformation, the closest nearest neighbor dis- 
tance is 2.63 A, which is larger than the 2.55 A value 
for TAO-2 and 2.57 A for Silcock. Because the nearest- 
neighbor distances in a and u> are 2.95 A and 2.65 A, 
respectively, it is not surprising that TAO-1 has the low- 
est barrier. The barrier to pass through the intermediate 
/3 phase, as suggested by Usikov and Zilbershtein 0, is 
108 meV/atom — drastically larger than TAO-1 's barrier 
of 9 meV/atom. 

Figure 0] shows the enthalpy along the a^ut transfor- 
mation path for the three mechanisms as a function of 
pressure. Our calculation shows u> slightly lower in en- 
ergy than a at GPa. The crystal structure at K has 
not been determined experimentally; however, extrapola- 
tion of the a-uj phase boundary indicates u) as the ground 
state 0. As pressure increases, the enthalpy of to rela- 
tive to a drops, and all three mechanisms decrease their 
enthalpy barrier. For the next four lowest landscape bar- 
rier mechanisms in Figure[21 we find GPa ab initio NEB 
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FIG. 4: Enthalpy barrier vs. pressure for the three lowest 
energy mechanisms, using ab initio NEB at 16 intermediate 
states. The ab initio NEB gives the T — energy along 
the homogeneous a—>ui pathway for each mechanism. The 
enthalpy barriers decrease with increasing pressure, (a) The 
TAO-1 mechanism gives the smallest barrier by a factor of 4 
at GPa. The TAO-2 mechanism (b) and Silcock mechanism 
(c) have larger barriers at GPa; an ordering which continues 
even at 40 GPa. 



barriers of 32, 37, 68, and 69 meV/atom. The barrier of 
the TAO-1 mechanism is lowest, even up to 40 GPa. 

We use classical nucleation theory to relate the homo- 
geneous barriers to the energy scales relevant for micro- 
scopic transformation mechanics, given by the formation 
energy of a critical nucleus. If the formation energy of 
the critical nuclei is of the same order as the barrier, our 
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mechanism prediction is not expected to change. In clas- 
sical nucleation theory, nuclei of new material grow only 
if they are larger than a critical size where the nucleation 
energy is an extrema. The nucleation energy consists of 
three terms: the interfacial energy, the enthalpy differ- 
ence, and the stress from the volume difference. 

In estimating the energy to create the critical nucleus, 
we use the a-Ti surface energy of 125 me v/A 2 as an up- 
per bound for the ot-uJ interfacial energy. The enthalpy 
difference at the experimental shock transition pressure 
of 10 GPa is —20 meV/atom, and the stress from the vol- 
ume difference is 3.5 meV/atom. We use the same crit- 
ical nucleus size for each mechanism in Table [I] because 
they each have a plane with no strain after the transfor- 
mation, so the nuclei likely grow as oblate spheroids to 
minimize stress. Classical nucleation theory estimates a 
critical nucleus radius of 240 A, (3.5 million atoms) with 
a formation energy of 30 keV. For that sized nucleus the 
9 meV/atom barrier of TAO-1 becomes a homogeneous 
barrier of 28 keV. In contrast, the competing TAO-2 and 
Silcock mechanisms have barriers of 200 and 120 keV re- 
spectively; far above the formation energy and the TAO-1 
barrier. 

Our systematic search, involving both tight-binding 
and ab initio calculations, combined with a nucleation 
analysis leads to the direct TAO-1 mechanism as the 
preferred a— >w mechanism in pure titanium. Our ap- 
proach allows us to find the lowest barrier mechanism for 
any martensitic phase transformation by generating and 
comparing all relevant possible mechanisms. Indeed this 
method provides an invaluable starting point for investi- 
gating any martensitic transformation. 
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